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Abstract. Many Galactic black hole and neutron star sources in low X-ray mass binaries show QPOs (quasi periodic oscilla- 
tions) in their observed X-ray fluxes, i.e., peaks in the Fourier variability power spectra. Pairs of twin peaks are observed, and 
in black-hole systems their frequencies Vupp, Vdown are in rational ratios. For example, in all four microquasars with twin peaks 
observed, v upp /v down = 3/2. The rational ratios have been postulated in a model that explained twin peak QPOs as a non-linear 
resonance between modes of acretion disk oscillations. For those microquasars where the mass of the X-ray source is known, 
we determine the black-hole spin, following from the observed QPO frequencies within various models of resonance. 
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1. Introduction 

Observed properties of the high frequency (fom ~1 Hz to ~10 3 Hz) quasi periodic oscillations (QPOs) in the X-ray variability of 
neutron star and black hole sources have been reviewed e.g. by van der Klis, (2000) and McClintock & Remillard (2003). There 
is no general agreement on a physical mechanism exciting QPOs. The resonance model proposed by Kluzniak & Abramowicz 
(2000) recalls these observed features that sharply illuminate the physical nature of the kHz twin peak QPOs: 

1. The Psaltis-Belloni-van der Klis-Mauche correlation, vi ow = 0.08vhi g h, between low and high frequency for all QPOs in 
black hole, neutron stars and white dwarf sources (see Figure proves that in general the QPOs phenomenon is due to 
accretion disk oscillations, and not to mere kinematic effects. It also proves that in general, QPOs frequencies do not scale 
with mass. However, QPOs in kilohertz range often come as twin peaks at frequencies v upp , Vdown- The twin peak kHz QPOs 
in microquasars seems to scale inversly with mass, v upp ~ l/M (see Figure This suggests a strong gravity origin of this 
particular class of accretion disk oscillations. 

2. In all four microquasars with twin peak kHz QPOs discovered, v upp /vd OW n = 3/2 (Table 1). This suggests that a non-linear 
resonance may be at work. 

With this in mind, it is natural to look for an explanation of twin peak kHz QPOs based on strong gravity 's oscillatons in a non- 
linear resonance. These are indeed two fundamental principles of the resonance model suggested by Kluzniak & Abramowicz. 

2. The strong gravity oscillations 

Consider a black hole 1 with the mass Mq and angular momentum Jq. In thin black hole accretion disks, matter spirals down the 
central black hole along stream lines that are located almost on the black hole equatorial plane 6 — 6q — tt/2, and that locally 
differ only slightly from a family of concentric circles r — ro — const. The small deviations, 5r = r-ro,50 — 6 — 6q are governed, 
with accuracy to linear terms, by 

8'r + co 2 6r = Sa r , 56 + a? e 58 = 6a e . (1) 

Send offprint requests to: G. Torok, terek@volny.cz 
1 We rescale mass with M = GM /c 2 = r G , and angular momentum with a = J c 3 /M 3 G 2 . We use Boyer-Lindquist coordinates, t, r, 8, 0. 
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Table 1. Frequencies of twin peak kHz QPOs in microquasars 



Source 


v upp [Hz] 


Vdown [Hz] 


2v U pp/3vdowi! _ 1 


Mass [ M e ] 


<a) XTE 1550-564 


276 


174 


0.05747 


8.4 — 10.8 


w GRO 1655^0 


450 


300 


0.00000 


6.0 — 6.6 


<a) GRS 1915+105 


168 


113 


0.00885 


10.0 — 18.0 


<b) H 1743-322 


240 


160 


0.00000 


not measured 



(a) From McClintock & Remillard (2003), (b) From Homan et al. (2003). 
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Fig. 1. (left) The vi ow = 0.08vhi g h correlation for neutron star binaries (open circles), black hole binaries (filled circles), and 
white dwarfs binaries SS Cyg (filled diamonds), and VW Hyi (open squares). Data for LMXRB from Belloni, Psaltis, & van der 
Klis (2002) who discovered the correlation for these sources, and for CV from Mauche (2002) who found that it also applies 
for CVs. See additional discussion by Psaltis, & van der Klis (2002) and Woudt & Warner (2002). Figure is reproduced from 
Mauche (2002) with permission of the Editors of the Astrophysical Journal. The correlation proves that the QPOs phenomenon 
must be due to accretion disks' oscillations and not to their kinematics like e.g. Doppler modulation of fluxes from isolated 'hot 
spots': QPOs are waves, not particles. Indeed, it is not possible that every time a spot is created at one radius to produce Vhi g h by 
Doppler modulation, another spot is created a t a precisely correlated larger radius to produce vi ow = 0.08vhi g h, also by Doppler 
modulation. Abramowicz, Barret, Kluzniak, Lee & Spiegel (2003) suggested that the low-high correlation in accretion disks 
could occur several reasons, well-known in oceanography. On deep waters, in a train of travelling waves with high frequency 
Vhigh, about every 9th wave is higher due to the side-band (Benjamin-Frei) instability. 9th waves cause a lower frequency feature 
vi ow that obviously obeys a correlation, vi ow = 0.11vhi g h- On shallow beaches, the lower frequency features in travelling waves 
are caused by backflows. Both effects could operate in accretion disks, producing the 1/0.08 » 13 th wave and the low-high 
correlation. 

Fig. 2. (right) The best fit for the frequency-mass scaling for microquasars' kHz twin peak QPOs, v upp = 2793/M [Hz/M ], 
adopted from McClintock & Remillard (2003) who pointed out the scaling. The scaling proves that oscillations are relativistic, 
i.e. they occur in strong gravity, at a radius fixed in terms of the gravitational radius rc = GM^/c 2 . 



Here dot denotes time derivative. For strictly Keplerian (i.e. free) motion, 8a r = 0, 6ag - 0, and the above equations describe two 
uncoupled harmonic oscillators with the eigenfrequencies cog, co r equal, in Kerr geometry (e.g. Nowak & Lehr, 1999), 



col = q k (l -4ax- y2 + 3a 2 X - 2 ), to 2 = Q| (l - 6 x" 1 + 8 a jT 3/2 - 3 a 2 x~ 2 ) , O k 



(GM \ m 



( x 3 ' 2 +a y 



(2) 



where x - r/M. The two epicyclic frequencies, vertical vg = cog/2n and radial v r = a> r /2n, are shown in Figure [3] together with 
the Keplerian orbital frequency vk = Q.kI^ for a nonrotating (a = 0) black hole, and for a moderately rotating (a = 0.8) black 
hole. Figures 0] and |5] show these frequecies in the whole range of a, from a maximally co-rotating black hole (a = 1) to a 
maximally counter- rotating black hole (a = -1). 
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Fig. 3. Orbital frequency v%, and the two epicyclic frequencies, radial v r , and vertical vg for Keplerian circular orbits around 
a 10 M black hole. Such orbits are posible only for radii larger than the radius of the circular photon orbit r p /,. This limit 
is labelled here and in other Figures by the subscript ph. Left panel for non-rotating black hole, right panel for a moderately 
(a = 0.8) rotating one. In Newton's theory with the 1/r potential, all the three frequencies are equal: vk — vg — v r . Strong 
Einstein's gravity makes vk > vg > v r . 

2.1. The highest possible orbital frequency: ISCO and RISCO 

In Newton's theory with the -GMo/r potential it is GMo/r 3 ' 2 = vk = v r = vg, but in the strong gravity of a rotating black hole, 
vk > vg > v r . The radial epicyclic frequency v r has a maximum at a particular circular orbit with the radius r > r ms , depending 
on the black hole spin (Kato & Okazaki, 1978), and goes to zero at r ms , the marginally stable circular orbit. The orbit is now most 
often called ISCO, which is a handy short name, and we suggest to call RISCO the marginally bound orbit at r m \,. 
Location of the inner edge of accretion disk, and therefore also the value of the highest possible orbital frequency, depends 
on disk's radiative efficiency. Very thin disks with high radiative efficiency have their inner edges almost exactly at ISCO. 
Correspondingly, the ISCO orbital frequency, vjc(/ mS ), is the highest possible orbital frequency for thin disks. However, as first 
proved analytically by Kozlowski, Jaroszyriski & Abramowicz (1978) and then confirmed in numerical simulations by many 
authors, for accretion disks with low radiative efficiency such as radiative and ion tori, slim disks, and adafs, the orbits are 
stabilized by the pressure gradient, and the inner edge goes close to RISCO. For such disks, the highest possible orbital frequency 
is obviously higher than the orbital frequency at ISCO, VK(r m b) > VK(/msX as illustrated in Figure 4 (right panel). 



2.2. Strong gravity's resonances and the inverse mass scaling 

Resonances occur when two of the three frequencies vk, v r , v$, are in rational ratios, e.g. v r /vg = n/m. Because all three frequen- 
cies vk, v r , vg have the general form, 



GMp 

, r „ 3 



Mr 



-1 



/(*, a) = 32293 I — I f(x, a) [Hz], (3) 



with f(x,a) being a dimensionless function, any particular resonance occurs at its own resonance radius x mm (a), and therefore 
each frequency at resonance scales approximately as 1/M. Radii of the three particular resonances (3:1, 2:1, 3:2) between the 
vertical and radial epicyclic frequencies are shown in Figure 6. 

3. The non-linear resonances 

3.1. Importance of the non-linear effects 

The effective potential H(r, 6; £) for orbital motion of a particle with a fixed angular momentum I > l ms has a minimum at ro({), 
corresponding to the location of a stable circular orbit. Its Talyor expansion (for simplicity we write it on the equatorial plane 
= tt/2), 

contains higher than quadratic terms, which means that small oscillations around the minimum at r - ro are described by non- 
linear differential equations (e.g. Landau & Lifshitz, 1976; Nayfeh & Mook, 1979). Non-linear resonances that may be excited 
in these non-linear oscillations have several characteristic properties that closely resemble those observed in QPOs: 
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Fig. 4. Left: The Keplerian orbital frequency for 1OM mass black hole. Right: RISCO and ISCO frequencies for 1OM mass 
black hole. The accretion disk inner edge must be located between RISCO and ICSO, depending on disk's efficiency. 
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Fig. 5. Left: The vertical epicyclic frequency for 1OM mass black hole. Right: the radial epicyclic frequency for 1OM mass 
black hole. 



1 . Resonance occurs in a region with a finite often large width 5v. 

2. The frequencies of oscillations v, depend on amplitude and for this reason they may be time dependent and may differ from 
the fixed eigenvalue frequencies v,(o) of the system, v,(f) = v,(0) + 5vi(t). 

3. Combination frequencies, v,-(o) ± V£(o) may be present. 

4. Subharmonic frequencies may be present. 



3.2. The 3:2 parametric resonance 

We shall start with an argument appealing to physical intuition and showing that the resonance now discussed is a very natural, 
indeed necessary, consequence of strong gravity. In thin disks, random fluctuations have 5r » 56. Thus, 5r56 is a first order term 
in 56 and should be included in the first order equation for vertical oscillations Q. The equation now takes the form, 

59 + u? g [1 + h 5r] 56 = 5a g , (5) 

where h is a known constant. The first order equation for 5r has the solution 5r = Aq cos(w r t). Inserting this in (0 together with 
5ag = 0, one arrives at the Mathieu equation (Ao is absorbed in h), 

59 + u? B [1 +h cos(Lo r t)]56 = 0, (6) 
that describes the parametric resonance. From the theory of the Mathieu equation one knowns that when 

— = — = -, n = 1, 2, 3..., (7) 

cog vg n 
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Fig. 6. Left: locations of the three resonances: the 3:2 parametric, and 2:1 and 3:1 forced for Schwarzschild case. Right: these 
locations depending on the black hole spin. Also shown is the radius rj? (dotted line) at which the standard relativistic Shakura- 
Sunyaev disk (Page & Thorne, 1974) emits locally the maximal flux, and the radius rc (dashed line) corresponding to the pressure 
centre of the maximally thick torus, i.e. the torus with constant angular momentum equal to the Keplerian value at RISCO. 

the parametric resonance is excited (Landau & Lifshitz, 1976). The resonance is strongests for the smallest possible value of n. 
Because near black holes v r < vg, the smallest possible value for resonance is n — 3, which means that 2 vg = 3 v r . This explains 
the observed 3:2 ratio, because, obviously, 

v upp = v g , v down = v r , (8) 

Of course in real disks neither 6r = Aq cos(co r t), nor 5ag = exactly, but one may expect that because these equations are 
approximately obeyed for thin disks, the parametric resonance will also be excited in realistic situations. And this is indeed 
the case. The parametric resonance of the type discussed above was found in numerical simulations of oscillations in a nearly 
Keplerian accretion disk by Abramowicz et al. (2003). Their numerical results were reproduced in an exact analytic solutions 
first by Rebusco (2003) and later confirmed and generalised by Horak (2003). The analytic solution is accurate up to third order 
terms in 6r, 66, and based on the method of multiple scales (see e.g. Nayfeh & Mook, 1976). Existence of the 3:2 parametric 
resonance is therefore a mathematical poroperty of thin, nearly Keplerian disks. It was found that the resonance is exited only 
in the non-Keplerian case, with some wea k forces Sag + and 6a r + present. Their origin is certainly connected to stresses 
(pressure, magnetic field, viscosity), but exact details remain to be determined — at present dag and 6a r are not calculated from 
first principles but described by an ansatz 2 . 

The parametric resonance occurs at a particular radius r^iia), determined by the condition 3co r {ry2, a) = Itogiryi, a) and equation 
(|2ji. We show the function r^ia) in Figure [6] In Figure we fit the 3:2 resonance theoretically predicted frequencies to the 
observational data for the three microquasars with the known masses. 



3.3. The forced 3:1 and 2:1 resonances 

A direct resonant forcing of vertical oscillations by the radial ones through a pressure coupling, and with Sa r ~ cos(a> r t), was 
evident in recent numerical simulations of oscillations of a perfect fluid torus (Lee, Abramowicz & Kluzniak, 2003). This supports 
an idea (Abramowicz & Kluzniak, 2001) for another possible model for the twin peak kHz QPOs: a forced non-linear oscillator, 

66 + co 2 g 68 + [non linear terms in 66] = h(r) cos(w r t), cog — n w r . (9) 

Obviously, there is no value for n such that cog and co r could be in the 3:2 ratio. However, the non-linear terms alow the presence 
of combination frequencies in a resonant solutions for 56( f) (see e.g. Landau & Lifshitz, 1974), 

CO- — a>g - OJ r , cj+ — cog + co r . (10) 

2 While the lack of a full physical understanding is obviously not satisfactory, the experience tells that such a situation is not uncommon for 
non-linear systems. Examples are known of mathematically possible resonances causing damage in bridges, areoplane wings etc., for which no 
specific physical excitation mechanism could have been pinned down (Nayfeh & Mook, 1976). 
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Fig. 7. Fit of the frequency v upp = vg(M, a) predicted by the 3:2 parametric resonance model to the frequencies observed in three 
microquasars with known masses. The spin parameter a is not known from observations, and the lines a = const are calculated 
from the model. One should notice that the deduced black hole spins are rather high, which is expected if microquasars' jets are 
a signature of large spin of the central black hole. 



Simple arthmetic shows (Abramowicz and Kluzniak, 2001) that these combination frequencies may be in the 3:2 ratio if and only 
if n = 2, or n = 3 and that in these two cases the observed frequencies Vd own = u)d 0wn /2n and v upp = a> up? /2n are uniquely given 
by, 

Wdown - to- - 2u> r , <jj upp = utg = 3 a> r for n — 3 forced epicyclic resonance tag = 3a> r , (11) 
w upp - to+ - 3u> r , Wdown = wg = 2u> r for n — 2 forced epicyclic resonance tog = 2w r . (12) 
We fit observed QPOs to these predicted by the forced epicyclic 3:1 and 2:1 resonances in Figure|S] 




3.4. Keplerian resonances 

The resonances discussed so far were due to a coupling between epicyclic oscillations — radial and vertical. Such a coupling 
exists in variety of realistic physical situations. It is more difficult to imagine a realistic physical coupling between radial epicyclic 
oscillations and orbital Keplerian oscillations. One possibility is connected to the fact that g-modes have pattern frequency 
2nv m = +(v r + mvjf), and these can be in co-rotation resonance, i.e., with v m = vk (Kato, 2001). The case VkI^t = 3/2 is 
excluded by observations (Figure 9, left panel). In addition, co-rotation resonance leads to damping, and not excitation, of modes 
(Kato, 2002). 

Another possiblity that may come to mind is based on the following idea (Spiegel, 2003). When the potential vorticity is con- 
served, coherent vortices tend to form in pairs with opposite vorticity (Bracco, Provenzale, Spiegel & Yecko, 1999). One may 
imagine that because the spatial distance between the two structures that oscillates with the epicyclic radial frequency, depends 
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on the velocity profile of the disk i.e. also on the oscillations of orbital velocity, a resonance between these two frequencies is 
possible. In Figure[9]we plot the prediction of the only three possible cases corresponding to w U pp/wdown = 3/2, 

Wupp = 3o> r , Wdown = Wk = 2(jJ r , 0) upp = 2 a>K, Wdown = <*>r = (1/3)«K, W upp = a>K, Wdown = <*>r = (2/3)«K- (13) 




4. Conclusions 

In Table 2 we summarize results of our estimates of the black hole spin for the three microquasars with known masses. The 
uncertainty in the spin estimate is due to uncertainty in the present knowledge of the black hole mass. In this respect we note that 
in the case of GRO 1655-40, the possible spin ranges given in Table 2 for 3:1 and 2: 1 Keplerian resonances are, respectively, 0.45 
- 0.53, and 0.31 - 0.42. They are more precise than the corresponding estimates by Abramowicz & Kluzniak (2001), i.e. 0.36 - 
0.67 and 0.2 - 0.6. The improvement in the accuracy is due to the improvement in the knowledge of the mass of GRO 1655-40. 
In 2001 the mass was known to be in the range 5.5 < M/M < 7.9, while in 2003 the accuracy increased to 6.0 < M/M Q < 6.6. 
The resonance model is rather sensitive to observational constraints: the data already excludes the 3:2 Keplerian resonance 
as a possible explanation of twin peak QPOs in case of two microquasars. All other resonannces discussed in the paper are 
consistent with the existing data, but it is plausible that future observations may narrow down the choice of a resonance. Future 
developements in accretion disk theory could also narrow down the choice. For example, it is often argued that the presence 
of relativistic jet is a signature of large black hole spin (Blandford and Znajek 1977, see however Fender et al. 2004). Because 
microquasars do have jets, one would except that a w 1 for their black holes. This argument, if proven true, would uniquely point 
to the 3:2 parameric resonance, as the only possible choice. 

Different resonances occur at very different resonance radii. Figure [6] shows that accretion disk physics at these radii is also 
very diferent. The 3:2 parametric resonance is located at a very outer part of the innermost region of the disk, the 2: 1 resonance 
is located in the middle of this region, and the 3:1 resonance most close to the inner edge. Therefore, a physical excitation 
mechanism (still unknown) must be very different for different resonances. 

If the resonace model is correct, and the kHz double peak QPOs are indeed due to non-linear strong gravity's resonance, the 
QPOs phenomenon would have a fundamental importance as a practical test of super-strong gravity. It could also be useful in 
several astrophysical applications of black hole astrophysics. For example, the l/M scaling of the twin peak QPOs frequencies 
with the 3:2 ratio (shown in Figure 1), was proposed by Abramowicz, Kluzniak, McClintock & Remillard (2003) as a method 
for estimating black hole masses in AGNs and ULXs, based on Mirabel & Rodriguez (1998) analogy between microquasars in 
our Galaxy and distant quasars. Indeed, if the analogy is also valid for accretion disk oscillations, then discovering in ULXs the 
twin peak QPOs frequencies with the 3:2 ratio, would resolve the controversy about their mass: if ULXs black holes have the 
same masses as microquasars, the frequencies will be ~ 100 Hz as in microquasars. If ULXs black holes are ~ 1000 times more 
massive, the frequencies will be ~ 0. 1 Hz instead. 
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